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Cell-cell and cell-matrix adhesions play essential roles in the function of tissues. There is grow- 
ing evidence for the importance of crosstalk between these two adhesion types, yet little is known 
about the impact of these interactions on the mechanical coupling of cells to the extracellular-matrix 
(ECM). Here, we combine experiment and theory to reveal how intercellular adhesions modulate 
forces transmitted to the ECM. In the absence of cadherin-based adhesions, primary mouse ker- 
atinocytes within a colony appear to act independently, with significant traction forces extending 
throughout the colony. In contrast, with strong cadherin-based adhesions, keratinocytes in a cohe- 
sive colony localize traction forces to the colony periphery. Through genetic or antibody-mediated 
loss of cadherin expression or function, we show that cadherin-based adhesions are essential for 
this mechanical cooperativity. A minimal physical model in which cell-cell adhesions modulate 
the physical cohesion between contractile cells is sufficient to recreate the spatial rearrangement of 
traction forces observed experimentally with varying strength of cadherin-based adhesions. This 
work defines the importance of cadherin-based cell-cell adhesions in coordinating mechanical activ- 
ity of epithelial cells and has implications for the mechanical regulation of epithelial tissues during 
development, homeostasis, and disease. 



Mechanical interactions of individual cells have a cru- 
cial role in the spatial organization of tissues [TJ [2] and 
in embryonic development [3H5]. The mechanical coop- 
eration of cells is evident in dynamic processes such as 
flow-induced alignment of vascular endothelial cells [6] 
and muscle contraction [7]. However, mechanical interac- 
tions of cells within a tissue also affect the tissue's static 
mechanical properties including elastic modulus [8 , sur- 
face tension [9], and fracture toughness [10] . Little is 
known about how these tissue-scale mechanical phenom- 
ena emerge from interactions at the molecular and cellu- 
lar levels [TT] . 

Tissue-scale mechanical phenomena are particularly 
important in developmental morphogenesis [12], home- 
ostasis [13], and wound healing [14] in epithelial tissues. 
Cells exert mechanical force on each other at sites of in- 
tercellular adhesion, typically through cadherins [T5| fT6] . 
as well as on the underlying extracellular matrix (ECM) 
through integrins [T7HT9] . Cadherin-based adhesions can 
alter physical aspects of cells such as the surface tension 
of cellular aggregates [20] and the spreading [21] and mi- 
gration [22] of single cells adherent to cadherin-patterned 
substrates. Integrity of intercellular adhesions may also 
contribute to metastatic potential [23 . We and others 



have shown that epithelial cell clusters with strong cell- 
cell adhesions exhibit coordinated mechanical behavior 
over length scales much larger than a single cell [24H26] . 
Several studies have implicated crosstalk between cell- 
ECM and cell-cell adhesions [27j [28] that can be modu- 
lated by actomyosin contractility [29]. Recent data sug- 
gest that integrin-mediated adhesions can modulate the 
composition [30j [31] and tension [25] [32j [33] of cell-cell 
junctions. While cadherins have been shown to modify 
local traction forces [34] and monolayer contractility [35] , 
the effects of intercellular adhesions on the spatial orga- 
nization of cell-ECM forces remain unexplored. 

In this paper, we address the impact of intercellular ad- 
hesions on cell-ECM traction forces in colonies of primary 
mouse keratinocytes. We measure tractions of colonies 
of keratinocytes before, during, and after formation of 
cadherin-mediated intercellular adhesions. As cadher in- 
dependent junctions form, there is dramatic rearrange- 
ment of cell-ECM traction forces from a disorganized, 
punctate distribution underneath the colony to an or- 
ganized concentration of force at the colony periphery. 
Through perturbations of cadherin-based adhesions, we 
demonstrate an essential role for cadherin in organizing 
cell-matrix mechanics. Finally, the spatial reorganization 
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of cell-matrix forces is reproduced by a minimal phys- 
ical model of a cell colony as two-dimensional objects 
connected by springs and adherent to a soft substrate. 
While downstream signaling pathways regulate responses 
to cadherin-based-j unction formation, our experimental 
data and physical model suggest that the simple physical 
cohesion created by intercellular adhesions is sufficient to 
organize traction forces. These results have implications 
for intercellular adhesions' role in the mechanical rela- 
tionship of tissues to their surroundings and the emer- 
gence of tissues' bulk material properties. 



RESULTS 

Traction Stresses Dynamically Reorganize in 
High-Calcium Medium 

To investigate the relationship between cadherin-based 
intercellular adhesions and cell-matrix traction stresses, 
we induced the formation of cadherin-based adhesions in 
primary mouse keratinocytes by elevating extracellular- 
calcium concentrations. In low-calcium medium, ker- 
atinocytes plated at low density proliferated into colonies 
of cells with weak cell-cell interactions. Exposing ker- 
atinocytes to high-calcium medium resulted in formation 
of cadherin-based cell-cell adhesions after 6-12 h. 

We quantified the effect of cell-cell adhesions on cell- 
matrix forces using traction force microscopy (TFM) [36 . 
We plated keratinocytes onto fibronectin-coated, elastic 
silicone gel coupled to glass. To quantify gel deforma- 
tion due to cell-ECM traction force, we imaged fluores- 
cent beads embedded in the silicone gel and measured 
the beads' displacements relative to their positions af- 
ter removing the cells with proteinase K. We calculated 
in-plane traction stresses, cr^, from bead displacements 
and the substrate's elastic properties [371 EH] (Materials 
& Methods). 

Over 12 h in high-calcium medium, keratinocytes de- 
veloped cell-cell junctions [39] and contracted [40] (Fig. 
[TJt4-(7). Prior to adhesion formation, in-plane traction 
stresses emanated from both the colony periphery and 
the interior junction of the three cells in a colony. Forces 
at the colony periphery pointed radially inward, while 
interior forces pointed in various directions (Fig. [l]D). 
During the timecourse, traction stress in the middle of 
the colony gradually weakened (Fig. [l]^), and by 12 h 
after calcium elevation, interior traction stress all but 
disappeared (Fig. [T]F). 

From substrate displacement and traction stresses, we 
calculated the strain energy density, w, the mechanical 
work per unit area performed by the colony to deform 
the substrate [41] (Materials & Methods). Shortly after 
calcium elevation, high strain energy was localized both 
underneath and at the periphery of the colony (Fig. [!]£?). 
12 h after calcium elevation, strain energy was limited to 
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FIG. 1: Traction stresses dynamically reorganize in high- 
calcium medium. (A—C) DIC images of a three-cell colony 
at 45 min (A), 6h (B), and 12 h (C) after calcium elevation. 
(D-F) Distribution of in-plane traction stresses (red arrows) 
for cell colony at timepoints in A—C overlaid on DIC im- 
ages. For clarity, one-quarter of calculated traction stresses 
are shown. (G-I) Distribution of strain energy density, w, for 
cell colony at timepoints in A—C. Blue lines mark individual 
cell boundaries. (J) Schematic for calculating azimuthal-like 
averages for strain energy. Colony outline is eroded inward 
by distance, A, in discrete steps, 5, until entire colony area 
has been covered. Average strain energy density is then cal- 
culated for each concentric, annular-like region. (K) Strain 
energy profiles for three-cell colony at six timepoints after 
calcium elevation. Solid colored lines represent colony's aver- 
age strain energy density as a function of distance, A, from 
colony edge. Each profile is mirrored about A « R, the ef- 
fective colony radius. Colony periphery (A = 0) is indicated 
by dashed vertical black lines. Strain energy at A < cor- 
responds to regions outside colony periphery. (L) Average 
strain energy density for entire colony at 15 min intervals from 
30 min to 12 h after calcium elevation. Plot colors in K and L 
are scaled according to time, t, after calcium elevation, from 
cyan at t = to magenta at t = 12 h. Scale bars in A-I 
represent 50 /xm. 



the colony edge (Fig. [Tjf). 

To quantify these spatial changes, we calculated 
azimuthal-like averages of strain energy during the time- 
course. We eroded the colony outline inward by distance, 
A, in discrete steps, 5, until the entire colony area was 
covered (Fig. [!}/). We calculated the average strain en- 
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ergy, w(A), in each of these concentric, annular-like re- 
gions and plotted it as a function of distance from the 
colony edge, A (Fig. [TJ/T). During the first 3h after cal- 
cium elevation, three peaks exist in the strain-energy pro- 
files, corresponding to localization of strong strain energy 
at the colony periphery (A = 0) and center. Between 5 
and 9h, the center strain-energy peak diminishes and 
disappears, and high strain energy is only at the colony 
periphery. We measured some strain energy outside the 
colony (A < 0) due to the finite spatial resolution of our 
implementation of TFM. 

Although strain-energy localization changed after cal- 
cium elevation, the colony's overall average strain energy 
density was relatively consistent during the timecourse 
(Fig. Hotspots of strong strain energy (yellow re- 

gions in Fig. [l]G) were no longer present by the end of 
the experiment (Fig. [l]^), but overall average strain en- 
ergy density was compensated by a decrease in colony 
area. 



Traction Stresses Systematically Reorganize in 
High-Calcium Medium 

To probe how intercellular adhesions alter traction 
forces across a large range of colony geometrical size and 
cell number, we analyzed the magnitude and localiza- 
tion of traction force in 32 keratinocyte colonies in low- 
calcium medium and 29 keratinocyte colonies after 24 h 
in high-calcium medium. A total of 117 low-calcium cells 
and 150 high-calcium cells comprised these colonies, each 
containing 2-27 cells, and spanned a geometrical dynamic 
range of nearly a factor of 100 in spread area. 

In general, low-calcium colonies exhibited traction 
stresses throughout the colony, usually pointing radially 
inward from the colony edge and in various directions 
in the interior (Fig. |2]A). Regions of high strain energy 
were found throughout the interior (Fig.[2]£?). In contrast, 
high-calcium colonies displayed traction stresses generi- 
cally pointing radially inward from the colony edge (Fig. 
[2]C) with hardly any strain energy beyond the colony 
edge (Fig. |D). 

To quantify these spatial distributions, we plotted av- 
erage strain energy density as a function of distance, A, 
from the colony edge (as depicted in Fig. [l|/). Average 
strain energy densities, w(A), were normalized by the av- 
erage strain energy density at the colony periphery, w(0). 
These profiles (Fig. [2] E and F) terminate where inward 
erosion covered the entire area of the colony, at A « i?, 
where R is the effective radius of the colony, given by the 
radius of the disk with the same area as the colony. 

In most low-calcium colonies, we observed some local- 
ization of strain energy at the colony periphery (A = 0) 
and high amounts of strain energy throughout the colony 
(A > 0), sometimes at the colony center (A « R) 
(Fig. In contrast, the strain energy of nearly all 



the high-calcium colonies was strongly localized at the 
colony periphery, generally decaying to zero toward the 
colony center (Fig. |5]F). Although this trend seems to 
hold regardless of number of cells in the colony, the dif- 
ference is much less pronounced for the smallest colonies 
(R < 50 /im). The radii of small colonies are compara- 
ble to the traction-stress penetration length, £ p , which 
measures how far from the periphery traction stresses 
penetrate the colony. Thus in small colonies, the stress 
measurements do not readily distinguish forces gener- 
ated at the colony center or periphery. In our previ- 
ous study on high-calcium keratinocytes, we measured 
£ p « 11 /im [26]. 

Next, we quantitatively compared the spatial distri- 
butions of strain energy across these two colony popula- 
tions with and without cadherin-based intercellular adhe- 
sions. We calculated the total strain energy, W exerted 
by each colony and the relative distance into the colony 
from its periphery, A/i?, required to capture 75% of the 
total strain energy, 3W/4. We separated larger colonies 
(R ^> £ p , or R > 50 /im) of the low- and high-calcium 
populations. Large, low-calcium colonies required on av- 
erage 10% more inward erosion (statistically significant, 
p = 0.0002) to achieve 75% of the total colony strain 
energy than large, high-calcium colonies, whereas there 
was no significant difference in strain energy distribution 
for the populations of small (R < 50 /im) colonies (Fig. 
|2]G) (p = 0.43). These data suggest that formation of 
cadherin-based adhesions in high-calcium medium results 
in a shift in localization of traction stress from internal 
regions of the colony to the periphery. 

The low- and high-calcium colonies did not seem to 
exhibit different amounts of average strain energy den- 
sity. A plot of total strain energy versus colony area, A, 
while scattered, shows no apparent difference between 
these populations (Fig. |2]/J). In general, larger colonies 
tended to perform more work on the substrate. 

Because low- and high-calcium keratinocyte colonies 
have different arrangements of cytoskeletal and adhesion 
proteins, we characterized spatial localizations of actin, 
E-cadherin-mediated cell-cell adhesions, and focal ad- 
hesions in keratinocyte colonies using phalloidin stain- 
ing and immunohistochemistry (Materials & Methods). 
E-cadherin is highly expressed in keratinocytes, mediates 
adhesive activity, and is essential for adherens-j unction 
formation. In high-calcium colonies, E-cadherin was lo- 
calized at keratinocyte junctions (Fig. [5]^). Positions of 
actin stress fibers were correlated with areas of strong 
E-cadherin localization, and there was coordination of 
the orientation of actin fibers across multiple cells, con- 
sistent with earlier reports on cytoskeletal rearrangement 
after calcium elevation [13]. While traction stresses of 
low- and high-calcium colonies had different spatial dis- 
tributions, focal adhesions, marked by paxillin, were con- 
centrated at the colony periphery in both cases. 
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FIG. 2: Traction stresses systematically reorganize in high- calcium medium. (A) Distribution of in-plane traction stresses 
(red arrows) of an eight-cell wildtype colony in low-calcium medium overlaid on DIC image of colony. For clarity, one-ninth 
of calculated traction stresses are shown. (B) Strain energy distribution, w, of low-calcium colony in C with individual cell 
outlines in blue. (C) Distribution of traction stresses (red arrows) of a six-cell wildtype colony in high-calcium medium for 24 h 
overlaid on DIC image of colony. For clarity, one-ninth of calculated traction stresses are shown. (D) Strain energy distribution, 
w, of high-calcium colony in E, with individual cell outlines in blue. (E) Strain energy profiles for n = 32 low-calcium colonies. 
(F) Strain energy profiles for n = 29 high-calcium colonies. In E and F, each solid curve represents a colony's average strain 
energy density as a function of distance, A, from colony edge. Each profile terminates where inward erosion covers entire colony 
area, at A « 7£, the effective colony radius, indicated by dashed line. The erosion is defined in legend of Fig. [Ij7. Average 
strain energy is normalized to value at colony periphery, w(0), giving each colony the same height on the graphs, indicated 
by the vertical scale bar. For clarity, profiles are spaced vertically according to colony size, with profiles for larger colonies 
(terminating at larger values of A) appearing higher up the y axis. Profile colors correspond to colony cell number given in 
the legend. (G) Quantification of relative distance from colony periphery (A/R) corresponding to 75% of total strain energy, 
3VF/4, in colonies in low- or high-calcium medium. Small colonies (R < 50 /mi, below hash marks in E and F), in low- (n = 8) 
or high-calcium (n = 8) medium showed no significant difference, whereas large (R > 50 /xm) low-calcium colonies (n = 24) 
had significantly more strain energy closer to colony center than large high-calcium colonies (n = 21). Statistical significance 
between low- and high-calcium populations is indicated by asterisks (p < 0.001). Error bars indicate one standard deviation. 
(H) Relationship between total strain energy, W, and area, A, of colonies in low- and high-calcium medium. Open symbols 
correspond to low-calcium colonies, closed symbols to high-calcium colonies. Symbol colors indicate colony cell number, given 
in the legend. (/ and J) Keratinocytes in low-calcium medium (7) or after 24 h in high-calcium medium (J) labeled with 
anti-E-cadherin and anti-paxillin antibodies and stained with phalloidin to mark F-actin. Scale bars in A—D, 7, and J represent 
50 /im. Data for high-calcium colonies in F-H are adapted from [26] . 



Cadherin-Based Adhesions are Required for 
Organization of Traction Stresses in High-Calcium 
Medium 

Because elevation of extracellular calcium modu- 
lates cellular properties in addition to cadherin-based- 



adhesion induction [42j [43], we sought to isolate the 
role of cadherin in spatially organizing traction forces. 
We used two different methods to inhibit formation of 
cadherin-based adhesions in the presence of high-calcium 




FIG. 3: Cadher in-based adhesions are required for organization of traction stresses in high-calcium medium. (A) Localization 
of E-cadherin, phalloidin (F-actin), and paxillin in colony of three wildtype keratinocytes in high- calcium medium for 24 h 
with DECMA-1. (B) Distribution of traction stresses (red arrows) of five-cell colony in high-calcium medium for 24 h with 
DECMA-1 overlaid on DIC image of colony. For clarity, one-sixteenth of calculated traction stresses are shown. (C) Strain 
energy of colony in B with individual cell outlines in blue. (D) Strain energy profiles for n = 15 DECMA-l-treated colonies. 
Each solid curve represents colony's average strain energy density as a function of distance, A, from colony the edge, as defined 
in Fig. [I]/. For clarity, profiles are spaced vertically according to colony size. Each profile terminates where inward erosion 
covers entire colony area, at A & R. (E) Localization of E-cadherin, phalloidin (F-actin), and paxillin in a colony of three 
E-cadherin-knockout/P-cadherin-knockdown (KO/KD) keratinocytes after 24 h in high-calcium medium. (F) Distribution of 
traction stresses (red arrows) of a colony of three KO/KD keratinocytes in high-calcium medium for 24 h overlaid on DIC 
image of colony. For clarity, one-sixteenth of calculated traction stresses are shown. (G) Strain energy distribution of colony 
in F with individual cell outlines in blue. (H) Strain energy profiles for n = 14 KO/KD colonies after 24 h in high-calcium 
medium. As in D, profiles were calculated as defined in Fig. [I]/. Profile colors in D and H correspond colony cell number given 
by the legend between Fig. ^E and F. (I) Comparison of the strain energy distribution for large low-calcium (n — 24), large 
high-calcium (n = 21), DECMA-1 (n = 15), and KO/KD (n — 14) colonies. Values represent proportion of the colony from 
periphery inward, A/R, necessary to capture 75% of total colony strain energy. Error bars indicate one standard deviation. 
Higher proportions indicate higher strain energy nearer colony center. Statistical significance between pairs of colony conditions 
is indicated by one asterisk (*) (p < 0.05), two asterisks (**) (p < 0.01), or three asterisks (***) (p < 0.001). Scale bars in A—C 
and E-G represent 50 /am. 



medium. First, we used the function-blocking antibody 
DECMA-1, which prevents homophilic binding between 
extracellular domains of E-cadherin [44]. DECMA-1 
was added to keratinocyte colonies with high-calcium 
medium for 24 h. Immunostaining of these colonies 
showed strong reduction of E-cadherin at intercellular 
contact (Fig.[3]A). Despite this change, we observed mini- 
mal coordination of actin across multiple cells in a colony, 
and focal adhesions were present at the colony periph- 
ery and throughout the colony interior. In keratinocytes 
in high-calcium medium with DECMA-1, we measured 
traction stress and strain energy throughout the colony, 
in particular at cell-cell contacts (Fig.[3]i? and C). Strain 
energy profiles of 15 DECMA-l-treated colonies (all with 
R > 50 //in) show many cases of high strain energy trans- 
mitted in the colony interior (Fig. |3]D). 

We further investigated the role of classical cad- 
herins using primary keratinocytes from an epidermal- 
E-cadherin-knockout (KO) mouse [45]. We used shRNA 
to knock down (KD) the other classical cadherin ex- 



pressed in these cells, P-cadherin, which is upregulated 
in E-cadherin- null cells [46] (Materials & Methods). We 
analyzed cell-cell and cell-matrix adhesions by immunos- 
taining KO/KD cells cultured in high-calcium medium 
for 24 h. Colonies of KO/KD cells showed no E-cadherin 
expression, did not coordinate their actin cytoskeletons 
across multiple cells, and displayed a slight reduction of 
focal adhesions underneath the colony (Fig.[3]E). As with 
DECMA-l-treated colonies, KO/KD colonies in high- 
calcium medium for 24 hr showed traction stresses and 
strain energy underneath cell-cell contacts (Fig.[3]F and 
G). Strain energy profiles of 14 KO/KD colonies in high- 
calcium medium (all with R > 50 jam) show strong strain 
energy transmitted throughout the colony (Fig. |3]/J). 

DECMA-l-treated colonies needed on average 6% 
more inward erosion than large high-calcium wildtype 
colonies to achieve 75% of the total colony strain energy 
(statistically significant, p = 0.048). KO/KD colonies 
required on average 10% more inward erosion than large 
high-calcium wildtype colonies to achieve 75% of the total 
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A [fj,m] A [fim] A [fim] 

FIG. 4: Minimal physical model captures cadherin- dependent 
organization of traction stresses. (A) Schematic of planar 
colony of three hexagonal cells. (B-D) Strain energy dis- 
tributions for colony of three hexagonal cells with different 
spring stiffness, k, expressed in units of E/L, where E is the 
Young's modulus of the cell and L the side length of each 
hexagon. (E-G) Spatial profiles of average strain energy as 
a function of distance, A, from colony edge for different val- 
ues of k corresponding to data in B-D. Other parameters: 
£ p /L 0.2, E lkPa, v 0.4, a a 4kPa, h 0.2 pm, 
Y — 2 x 10 6 N/m 3 {Appendix). Scale bars in B-D represent 
50 /im. 

colony strain energy (statistically significant, p = 0.002) 
(Fig. [3]^). Compared to large low-calcium colonies using 
this same measure, neither DECM A-l-treated colonies 
(p = 0.14) nor KO/KD colonies (p = 0.94) showed signif- 
icant differences in spatial distributions of strain energy. 
Thus, keratinocytes in high-calcium medium are unable 
to organize traction forces to the colony periphery in the 
absence of cadherin-based cell-cell junctions. 

Minimal Physical Model Captures 
Cadherin-Dependent Organization of Traction 
Stresses 

Because of the simple spatial trends of traction stresses 
observed in colonies with and without intercellular ad- 
hesions, we examined whether a minimal physical model 
could reproduce the experimental results. We model each 
cell in a colony as a homogeneous and isotropic elas- 
tic material [47J 08]. In our model, each cell exerts a 
contractile "pressure" opposed by strong adhesion to a 
compliant substrate [49]. At each point within a cell, we 
require that these opposing forces balance. This model 
ignores all active processes modulated by cell-cell ad- 
hesions, including downstream signaling, and represents 
each intercellular adhesion as a purely physical connec- 
tion characterized by a spring constant, k [50] . 

To make predictions with this model, we use a numer- 
ical solution of the two-dimensional governing equations 
(Appendix). To mimic the cell geometry in the time- 
course experiment (Fig. [T]), we consider the case of three 
hexagonal cells (Fig. |4]A). We find that, for increasing 



cell-cell-coupling strength, fc, traction stress and strain 
energy disappear under cell-cell junctions (Fig. [I] B-D) , 
recapitulating the transition seen in real cells stimulated 
by calcium elevation (Fig. [I] D-F). The similarity be- 
tween model and experiment is also evident in plots of 
strain energy density as a function of distance from the 
colony edge (Fig.g^-G and Fig.^). 

The model demonstrates the importance of 
intercellular-adhesion strength in spatially organiz- 
ing cell-ECM forces. For weak cell-cell coupling (small 
fc), individual cells deform the substrate independently of 
each other, with significant substrate deformation at all 
edges of each cell. On the other hand, strongly coupled 
colonies (large k) behave as a cohesive, contractile unit, 
with substrate deformation only at the colony periphery. 

This planar model is an extension of an analytically 
tractable, one-dimensional model (Appendix). 

DISCUSSION 

Our results show that cadherin-based cell-cell adhe- 
sions modulate force transmission to the ECM. In partic- 
ular, our tract ion- force data on cohesive cell colonies sug- 
gest that intercellular-adhesion formation through classi- 
cal cadherins reorganize the spatial distributions of trac- 
tion stress. In colonies of cells with strong E-cadherin- 
based adhesions, cell-ECM traction stresses are localized 
in a ring around the colony periphery. In weakly cohesive 
colonies, regions of high traction stress appear through- 
out the colony. Furthermore, traction stresses cannot re- 
organize in high-calcium medium when cadherin-based 
adhesion is inhibited. Comparison of our experimen- 
tal data with our minimal physical model suggests that 
strong physical cohesion between cells is sufficient to 
drive the relocalization of cell-ECM forces to the periph- 
ery of cell colonies. While our data show that E-cadherin 
is necessary to reorganize traction forces, E-cadherin 
alone may not be sufficient. Further study is required to 
determine whether additional adhesive processes down- 
stream of adherens junctions, such as the formation of 
desmosomes by nonclassical cadherins [46 , are necessary 
to achieve sufficient cohesion. 

Our findings resonate with recent studies on cellu- 
lar adhesion pointing toward crosstalk of cadherin- and 
integrin-based adhesions. Focal adhesions have been ob- 
served to disappear underneath cell-cell contacts [30j[5T], 
but this effect may depend on substrate stiffness [28] 
and the extent of cell spreading [52]. Recent work has 
also suggested that forces transmitted through focal ad- 
hesions can modulate intercellular forces [25 , 28J, which 
in turn can modulate intercellular-junction assembly and 
disassembly [15] . Our study highlights intercellular adhe- 
sions' ability to impact cell-ECM force generation, which 
allows for bidirectional feedback between cell-cell and 
cell-matrix forces. Indeed, tension at cadherin junc- 
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tions p~3j [53] is known to elicit cell-signaling events and 
actin dynamics [5TJ [54H57] and contribute to collective 
cell migration [581 EH]. In light of these prior results 
on integrin-cadherin feedback, it is somewhat surpris- 
ing that a minimal physical model can capture the ob- 
served dependence of cell-matrix forces on the strength 
of cadherin- mediated cell-cell adhesions. 

Reorganization of cell-ECM forces is likely one impor- 
tant mechanism by which cadherin-based adhesions drive 
tissue morphogenesis and homeostasis. In development, 
differential adhesion has been shown to play an impor- 
tant role in cell sorting [60H62] , and the reorganization of 
intercellular forces in this context is entirely unexplored. 
Furthermore, in wound healing, we expect strong cell- 
ECM forces to be generated at a wound edge due to the 
local loss of intercellular adhesion. These forces could 
act as a signal, inducing migratory behavior in epithelial 
cells [241 [63] , activating responses of stromal cells, and 
organizing the ECM [64-66 . A key avenue for future in- 
vestigations will be to explore how organization of force 
stimulates cellular responses within tissues. 
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MATERIALS & METHODS 



Preparation of Substrates for Traction Force Microscopy. A 

borate buffer solution was made from deionized water with 3.8mg/ml 
sodium tetraborate and 5mg/ml boric acid. Silane (3-aminopropyl 
triethoxysilane) (Polysciences) was vapor-deposited onto 35 mm glass- 
bottom dishes (WillCo Wells) to allow florescent beads to be bonded 
to the surface. Beads were deposited by filling the dish with a so- 
lution containing dark- red florescent (660/680) carboxylate-modified 
microspheres with radius 0.1 /im (Life Technologies) at a volume ratio 
of 1:3,000 and 1 wt% l-ethyl-3-(3-dimethylaminopropyl)carbodiimide 
(EDC) (Sigma- Aldrich) at a volume ratio of 1:100 in borate buffer. 
Silicone elastomer was then prepared by mixing a 1:1 weight ratio of 
CY52-276A and CY52-276B (Dow Corning Toray). After being degased 
for 10 min, the elastomer was spin-coated onto the glass of the dish at 
2,000 rpm for 60s. The dish was cured at 50° C for 3 min and resulted 
in an elastic film ~21 /im thick. With the elastomer cross-linked, silane 
was vapor-deposited on the elastomer-coated dish. A second layer of 
florescent polystyrene beads was deposited at a higher concentration, 
volume ratios of 1:1,000 beads and 1:100 EDC in borate buffer. A sec- 
ond layer of fresh, degased elastomer was spin-coated at 10,000 rpm 
for 90 s resulting in a layer ~3 /im thick. The sample was cured at 
RT overnight. We estimated the Young's modulus, E, of the cured 
elastomer to be ~3kPa using bulk rheology. Before cells were plated, 
the elastomer surface was coated with fibronectin from bovine plasma 
(Sigma-Aldrich) at a concentration of 0.2mg/ml, which sat for 20 min 
at RT before being washed off with PBS. 

Confocal Microscopy. Images for TFM experiments were acquired 
using an Andor Revolution spinning-disk confocal system (Andor Tech- 
nology) mounted on an inverted microscope (Nikon Eclipse Ti) with a 
Plan Apo 60 X water-immersion objective lens with numerical aperture 
of 1.2 (Nikon). A 640 nm laser and DIC channel were used to image 
florescent beads and cells, respectively. Images were acquired with an 
iXon EMCCD camera with a resolution of 1,024 x 1,024 pixels (Andor 
Technology). The field of view was 113 X 113 fim 2 . Because a single 



field of view was too small to image an entire cell colony, between 9 and 
42 fields of view per colony were acquired, with adjacent fields of view 
overlapping by ~25% and stitched together with sub-pixel precision 
by aligning bead positions in overlapping regions. The stage was con- 
trolled through Motion Controller/Driver SMC100CC high-speed mo- 
torized actuators (Newport). We imaged fluorescent beads with con- 
focal image stacks of total thickness 5 /im to cover the beads' entire 
point-spread function in z. Confocal image slices were spaced 200 nm 
apart. These stacks were reduced to single images for particle-tracking 
by averaging the slices from five below to five above the slice with the 
highest total intensity. 

Live-Cell Imaging. Confocal image stacks of the fluorescent beads 
were acquired for each cell condition. Control images of the beads 
in their unstressed state were acquired after removing the cells from 
the elastomer with proteinase K (Life Technologies) at 0.5mg/ml for 
5 min and then washing with PBS. The cells on the microscope were 
maintained at 37° C using a heated microscope stage. pH was con- 
trolled using HEPES solution at 15 mM (Sigma-Aldrich). To inhibit 
the formation of cadherin-based adhesions, we added anti-E-cadherin 
antibody DECMA-1 (Abeam) at 6/2g/ml to the high-calcium medium. 
For consistency across cellular conditions, we controlled for colonies 
that deviated significantly from disk-shaped and contained cells with 
long protrusions by selecting for colonies whose actual perimeter, P, 
was no more than 1.5 times the perimeter of a circle of the same area, 
A, as the colony (P < 3V^A). 

Calculation of Traction Stresses and Strain Energies. Af- 
ter determining bead positions using centroid analysis in MATLAB |67| , 
we calculated the deformation of the substrate, u^(r,z ), across its 
stressed (with cells) and unstressed (with cells removed) states, where 
z is the distance between the substrate bottom and the bead layer. In 
Fourier space, the deformation field is related to the traction stresses 
at the surface of the substrate, h s , via linear elasticity, cr| z (k,/i s ) = 
Qij (k, z Q , h s )u s - (k, z Q ), where k represents the in-plane wave vector. 
Here, <r^ z (k, h s ) and n^(k,2 ) are the Fourier transforms of the trac- 
tion stress on the top surface and the displacements of the bead layer 
just below the surface, respectively. The tensor, Q, depends on the 
thickness and modulus of the substrate, the location of the beads, 
and the wave vector 37 38] . We calculated the strain energy density, 
u>(r) = \<J s iz (r, h s )u^ (r, h s ) 41 . The deformation on the surface was 
determined using u|(k, h s ) = QT. 1 (k, h s , h s )Q jk (k, z , h s )uf, (k, z ). 
Because of their small size and immersion in a viscous medium, we 
expect the colonies of cells to be in mechanical equilibrium (net force 
of zero). Due to experimental error in determining substrate displace- 
ment fields, we occasionally calculated non-zero net traction forces on 
a colony. We discarded colonies with more than 15% residual force, 

\f dA + ° S yz y) \ > 0.15/cL4|<7* z x + <7* z y|. 

Primary Keratinocyte Culture. Primary wildtype keratinocytes 
were isolated as described 68 . Briefly, isolated backskin of newborn 
CD1 mice was floated on dispase overnight at 4°C. The epidermis 
was separated from the dermis with forceps and incubated in 0.25% 
trypsin for 10 min at RT. Individual cells were released by tritura- 
tion and plated on mitomycin-C-treated J2 fibroblasts in low-calcium 
medium (0.05 mM CaCl2). After 2—4 passages, cells were plated on 
plastic dishes without feeder cells. Primary keratinocytes were also 
isolated as described |69| from newborn epidermis in which E-cadherin 
was conditionally deleted as described |45| . KO/KD cells were gen- 
erated by lentiviral transduction of E-cadherin-deficient keratinocytes 
using shRNA directed against P-cadherin, as described 46 . Cadherin- 
junction formation was induced by raising the concentration of CaCl2 
of the low-calcium medium to 1.5 mM. 

Immunohistochemistry. Cells were fixed in 3.7% formaldehyde 
for 10 min and then washed twice for 2 min in PBS. A blocking solution 
of normal goat serum, normal donkey serum, bovine serum albumin, 
gelatin, and triton X in PBS was used to prevent non-specific binding. 
Cells were stained using 3,000 units//il Alexa Fluor 594 phalloidin (Life 
Technologies) and primary antibodies 8ng//xl monoclonal mouse anti- 
E-cadherin (TaKaRa) and 4ng//xl rabbit anti-paxillin (Sigma-Aldrich). 
After being washed in PBS, cells were incubated with secondary anti- 
bodies 8ng/yiil goat anti-rabbit Alexa Fluor 488 (Life Technologies) 
and 8ng//xl goat anti-rat Alexa Fluor 647 (Life Technologies) and 
again with 3,000 units//il Alexa Fluor 594 phalloidin. Cells were then 
mounted in ProLong Gold with DAPI (Life Technologies). 

Fluorescent images of immunohistochemical staining were acquired 
using confocal laser scanning microscopy on a Zeiss LSM 510 system 
equipped with Ar, HeNe 543, and HeNe 633 laser lines allowing imaging 
with lasers of wavelengths 488, 568, and 633 nm and a Plan Apo 40 X 
oil-immersion objective with numerical aperture 1.3 (Zeiss). The field 
of view was 313 X 313 /im 2 with a maximum resolution of 2,048 X 2,048 
pixels. The stage was controlled using an MCU 28 unit (Zeiss). 

Statistical Analyses. Statistical significance for strain-energy dis- 
tributions was assessed with p values determined by two-sided Student's 
t tests. Statistically significant p values were those lower than 0.05. 
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FIG. 5: Minimal one-dimensional picture of N cells adhering 
via cadherin-based adhesions, modeled as linear springs of 
stiffness k. 



APPENDIX 

Planar Model of Cell Colonies as Elastic Media 

We model each cell in a colony as a contractile elastic ob- 
ject, assumed to be homogeneous and isotropic with Young's 
modulus, E, and Poisson's ratio, v |47l|48]. The constitutive 
relation for the thickness- averaged cellular stress tensor, cr^, 
with z, j denoting in-plane coordinates, is given by 

o-ij = 2 ^ E + ^ V ' u 6i i + diU 3 + d i u ^j + aaSi i' ( 1 ) 

where u is the cellular displacement field, Sij is the Kronecker 
delta, and o~ a > denotes the active stress due to actomyosin 
contractility. In the thin- film limit, force-balance requires 

hdjGij — Ym, (2) 

where h is the average cell height, and Y is a parameter de- 
scribing the strength of cell-substrate coupling and depends 
on substrate stiffness as well as on the strength of focal adhe- 
sions [T7lH9ir70] . 

We model cell-cell interactions as linear springs [50] of 
spring constant k per unit area, exerting a harmonic force 
f per unit area normal to the interface between two cells. 
The addition of springs translates into boundary conditions 
at the intercellular interfaces as crijUj — fi, with n denoting 
the outward unit normal. The edge of the colony, however, 
respects the stress-free boundary condition, aijUj = 0. We 
numerically solve Eqs. |l]) and |2]) subject to the aforemen- 
tioned boundary conditions using the matlab PDE Toolbox. 
We evaluate strain energy density, w, given by w = |T • u, 
where T = Yu is the local traction stress exerted by the 
colony. According to this model, spatial variation of trac- 
tion stresses and strain energy densities in the direction nor- 
mal to the edges is cha racterized by a penetration length, 
£ p = ^Ehv/[Y{l-v*)]. 



Model of Cell Colonies as Elastic Media in One 
Dimension 

As in the planar case, individual cells are described in one 
dimension as thin active elastic materials adherent to an elas- 
tic substrate. We consider N cells, each of rest length L/N 
and average height h, with cell-cell adhesions modeled by lin- 
ear springs of stiffness k (Fig. [5]). Let denote the internal 
stress in the a th cell and the corresponding displacement 
field. The one-dimensional equivalents of Eqs. |l]) and ^ for 
the a th cell are given by 

a (a) (x) = Bd x u (a) +a a (3) 

and 

hd x a (a) =Yu (a \ (4) 

respectively, where B is the longitudinal elastic modulus of 
the cell. In this one-dimensional picture, sub strate- induced 
traction penetration length is given by t v — \J Bh/Y. 

Internal stress distribution in the colony is then governed 
by equations 



f p dla (a) (x) + G a = a (a) (x) for 1 < a < N (5) 
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FIG. 6: (A) Internal stress and (B) strain energy density in a one-dimensional adherent cell-pair for kL/a a = 0.004 (dotted), 
kL/cfa = 0.4 (dashed), and kL/a a = 40 (solid). Parameters: £ P /L = 0.2, B/a a = 2, h/L = 0.1. (C) Intercellular force, / 
versus intercellular adhesion strength, k (in units of cr a /L) for £ P /L — 0.2. Inset: / as a function of £ P /L for kL/a a — 10. 



subject to boundary conditions 



<7 (1) U =0 =0, 

(a) I (a + l)| 

<? \x = cxL/N — CT \x = uL/N 



) (a) 



(6a) 



k 



y.L/N 



for 1 < a < N, 



cr y J \x=L = 0. 



(6b) 
(6c) 



For simplicity, we assume that the cell-cell adhesion springs 
have zero rest length and that the colony ends (x — 0, L) 
respect stress- free boundary conditions. 

Explicit solutions for cellular stresses in an adherent cell- 
pair (N = 2) are given by 



a^\x)/a a = 1 - exp 



(7a) 



2 sinh 



2k£ T 
B 



jt) [^ + exp(^)-l] ; 
[l + exp(^)]+exp(^)-l' 



cr^ 2 \x)/cFa — 2 sinh 



L — x 
2£r> 



(7b) 



cosh 



(^)-cosh(^f) + ^smh(^) 
^ cosh (^)+ sinh (± 

For weak intercellular coupling, k <C B/£ p , internal stresses 
are maximal at the center of individual cells and negligible at 
the cell-cell junction. For a strongly coupled cell-pair, k ^> 
B/£ p , internal stresses build up at the junction between the 
cells, which corresponds to the limit of a cohesive cell colony 
(Fig. 1614). In this case, internal stress takes the simple form 



a(x) 



cosh 



( L-2x \ 

V 2 *p J 



cosh 



(8) 



Strain energy density, w, in the cell-pair is determined using 
w(x) = ^T(x)u(x), where the traction T(x) = Yu(x). For a 
weakly coupled cell-pair (k — > 0), w is localized at the edges 
of individual cells, and the net traction force on each indi- 
vidual cell vanishes. In contrast, a strongly cohesive cell-pair 
(k — >• oo) behaves as a single cell, with strain energy density 
localized at the edge of the pair and vanishing at the junction 
(Fig. [6]S). For intermediate strengths of cell-cell adhesion, 
there is finite but small strain energy at the junction com- 
pared to the edges of the cell-pair. Traction force imbalance 
at each cell gives the estimate of the total force, /, transmit- 
ted to the intercellular adhesion, 



fL/2 



dx T (1) (x) 



h\a(0)-a(L/2)\ 

ha a |l - cosh (2^) I 



cosh 



(4>) 



(9) 



" haa k + B/2i/ ° T L>>£ *>- 

Intracellular force, /, grows monotonically with adhesion 
strength, k, before reaching a plateau when k is large 
(Fig. [6]C). However, the dependence of / on penetration 
length, £ p , which is inversely related to substrate stiffness, 
is non-monotonic (Fig. |6]C, inset). This biphasic relation 
arises from the competition among different elastic compo- 
nents (cell, substrate, and the intercellular spring) connected 
in series. For small £ p , the substrate is deformed less com- 
pared to the cells and to the intercellular spring, leading to 
a rise in intercellular force. A more compliant substrate with 
large £ p is likely to accommodate larger cellular forces, reduc- 
ing the net force transmitted to the intercellular adhesion. 



